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ABSTRACT 

Using the 100-m Effelsberg radio telescope operating at 1.36 GHz, we have per- 
formed a targeted radio pulsar survey of 289 unassociated 7-ray sources discovered by 
the Large Area Telescope (LAT) aboard the Fermi satellite and published in the 1FGL 



catalogue ( Abdo et al. 2010a]). This survey resulted in the discovery of millisecond pulsar 



J 1745+ 1017, which resides in a short-period binary system with a low-mass companion, 
M Cjm i n ~ 0.0137Mq, indicative of "Black Widow" type systems. A two-year timing 
campaign has produced a refined radio ephemeris, accurate enough to allow for phase- 
folding of the LAT photons, resulting in the detection of a dual-peaked 7-ray light-curve, 
proving that PSR J1745+1017 is the source responsible for the 7-ray emission seen in 
1FGL J1745.5+1018 (2FGL J1745.6+1015; |Nolan et al.|20T2] >. We find the 7-ray spec- 



trum of PSR J1745+1017 to be well modelled by an exponentially-cut-off power law 
with cut-off energy 3.2 GeV and photon index 1.6. The observed sources are known to 
contain a further 10 newly discovered pulsars which were undetected in this survey. Our 
radio observations of these sources are discussed and in all cases limiting flux densi- 
ties are calculated. The reasons behind the seemingly low yield of discoveries are also 
discussed. 

Key words: pulsars: general - pulsars: individual: PSR J1745+1017 - gamma-rays: 
general 
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1 INTRODUCTION 

The detection of pulsed 7-ray emission from the Crab Pulsar in 
the early 1970's ( TVasseur et al.|1970[|Grindfay|1972"[ >, the first 
of its kind, brought new light to the study of pulsar emission 
physics and high-energy emission physics in general. Gamma- 
ray photons of energies greater than 100 keV are created in 
processes involving nuclear or other non-thermal reactions, and 
as such become important when exploring the Universe at its 
most energetic. The current model for the creation of 7-ray pho- 
tons that we see from pulsars is that charged particles stripped 
from the surface are accelerated to relativistic energies in the 
pulsar's strong electric field. As these particles travel along the 
curved magnetic field lines, they produce 7-ray photons via syn- 
chrotron radiation, curvature radiation (e.g. Ruderman & Suther- 
land 1975 ) and inverse Compton scattering from lower-energy 
photons (e.g. |Daugherty & Harding 1986). The study of these 
processes gives insight into the structure and composition of the 
magnetospheres of pulsars. 

Prior to 2008, the most successful space-based 7-ray ex- 
periment was the Compton Gamma-Ray Observatory (CGRO), 
which was in orbit for nine years and carried the Energetic 
Gamma-Ray Experiment Telescope (EGRET, Kan bach et al.| 
1989). EGRET was sensitive to 7-ray photons in the range 20 
MeV - 30 GeV, and during its lifetime brought the known num- 
ber of 7-ray emitting pulsars up to at least six ( Thompson 2008 1. 
However, the legacy of EGRET for the radio community was not 
the pulsars it detected, but rather those sources for which it could 
make no positive association. Targeted radio searches of these 
169 7-ray sources, unassociated with either pulsars or blazars, 
were performed, leading to several pulsar discoveries (e.g. Keith 
|et al.|2008[ [Champion et al.|2005) . 

The Large Area Telescope (LAT) JAtwood et "aL][2009| ) 
aboard the Fermi Gamma-ray Space Telescope, represents a 
significant improvement upon EGRET, providing a greater en- 
ergy range and sensitivity, allowing for better measurements 
of source characteristics and localisations. With a host of 
new sources discovered, including many active galactic nuclei 
(AGNs) and pulsars, the Fermi LAT is the most successful GeV 
7-ray observatory to date. As with EGRET, it is those sources 
for which Fermi cannot immediately provide an association that 
have piqued the interest of the pulsar searching community. 
A catalogue of 1451 7-ray sources detected above 100 MeV 
was created from the first 11 months of LAT data. Of these 
sources, 630 were unassociated with known astrophysical ob- 
jects (AGNs, pulsars, etc.; Abdo et al. 2010a). Multi-wavelength 
observations of the unassociated sources were encouraged so as 
to determine their natures, with many radio observatories search- 
ing for radio pulsations in the Fermi observational error ellipses 
(e.g. |Keith et al.|201 1 [[Ransom et al.|201 l||Cognard et al.|201 1| >. 

While Fermi LAT data have already been proved to contain 
a wealth of pulsars, with more than 100 pulsars detected through 
blind periodicity searches and phase folding of LAT photons us- 
ing known pulsar ephemerides ( |Ray & Saz Parkinson|[20TT| ), 
low photon counts introduce strong selection biases in the de- 
tection of pulsars through blind searches of the LAT data. This 
is due to the large amount of computation required to perform 
wide-parameter-space searches of sparse photon data sets. For 
this reason, blind searches of the LAT data currently have great 



difficulty in detecting millisecond pulsars (MSPs) or pulsars in 
binary systems. 

Radio pulsation searches are subject to different biases and 
thus are an important alternate method for identifying LAT unas- 
sociated sources as pulsars. At the time of this writing, there 
have been 47 radio-loud pulsars discovered through searches of 
these sources, of which 41 are MSPs likely to be associated with 
their corresponding LAT source ( |Ray et aT||2012[ >. These dis- 
coveries highlight the importance of targeted radio searches of 
the LAT data, as these pulsars were most likely undetected in 
more general surveys due to shorter integration times or lack of 
searching for binary motion. 

Of the MSPs discovered, 10 are thought to be in 'Black 
Widow' systems where the companion star has a very low mass 
due to ablation by the strong wind of the pulsar ( |Fruchter et al.| 
1988). Before Fermi only three of these systems were known 
to exist outside globular clusters (Fruchter et al. 1988; Stappers 
et al. 1996; Burgay et al. 2006), which stresses the importance 
of investigating this new population of pulsars uncovered by the 
LAT. Those MSPs discovered which are not in Black Widow 
systems may also be of great use to current and future pulsar 
timing arrays for gravitational wave detection (Fost er & Backerj 
1 990 1, which benefit from an even distribution of precisely timed 
pulsars across the sky. 

In this paper we present a targeted search of 289 unassoci- 
ated Fermi sources using the 100-m Effelsberg telescope oper- 
ating at 1.36 GHz. The search has resulted in the discovery of 
a 2.65-ms pulsar, PSR J1745+1017, in a 17.5-hour binary orbit 
with a 0.016-Mq companion. The positions of 10 pulsars found 
in other targeted searches of unassociated LAT sources are con- 
tained within the 289 sources observed. For these sources, we 
discuss possible reasons for our non-detections and provide flux 
density limits where applicable. 

This paper is structured as follows. In Section[2] we discuss 
selection criteria for sources to be observed. In Section[3]we dis- 
cuss the observational methods and data processing. In Section 
[4]we discuss the survey sensitivity. In Section[5]we discuss sim- 
ulations of the survey. In Section [6] we discuss the results of the 
survey. In Section [7] we discuss the source selection and detec- 
tion rate. In Section[S]we present our conclusions. 



2 FERMI CATALOGUE SOURCE SELECTION 

All the sources searched for this paper were selected either 
from the Fermi LAT First Source catalogue (1FGL) (Abd o et al.| 
2010a ) or from an unpublished update to the 1FGL which cov- 
ered the first 18 months of Fermi observations. 

The 1FGL contains 630 sources unassociated with any 
known astrophysical object, in the 100 MeV to 100 GeV en- 
ergy range. The catalogue includes source localisations defined 
in terms of an elliptical fit to the 95% confidence level, power- 
law spectral fits, monthly light curves and flux measurements in 
5 energy bands for each source. A full description of the 1FGL 
can be found in Abdo et al. ( 2010a). To narrow down the num- 
ber of potentially observable sources, those below —20° decli- 
nation, and those for which a known association existed were 
ignored. For all sources, the pointing position was chosen to co- 
incide with the centre of the error ellipse. As better source local- 
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isations became available, pointing positions were altered such 
that they coincided with the center of the updated error ellipse. 

The remaining sources were ranked using an application 
of the Gaussian-mixture model as outlined in |Lee et al. | ( |2012^ . 
The best pulsar candidates exhibit significantly curved emission 
spectra and little 7-ray flux variability over time, in contrast with 
e.g. blazars and other AGNs (see e.g. Figure 17 of |Nolan et al.| 
[20T2> . 

The top 250 ranked candidates from the 1FGL catalogue 
were selected to be observed. A further 39 highly ranked sources 
from the 18-month update to the 1FGL were also selected to be 
observed, giving a total of 289 target sources. 



3 OBSERVATIONAL METHOD AND DATA 
PROCESSING 

All search data presented in this paper were taken with the 100- 
m Effelsberg radio telescope at a centre frequency of 1.36 GHz 
between November 2009 and July 2010, using the central horn 
of the new Effelsberg multi-beam receiver with 300-MHz band- 
width. 

The 289 sources selected for this survey were observed in 
three different observing campaigns. Initially all sources were 
observed with between 10- and 16-minute integrations. These 
integrations allowed for a preliminary shallow sweep of all 
sources such that the brightest pulsars, those likely to have been 
of the greatest use for timing applications, could be found. The 
second observing campaign was comprised of 32-minute inte- 
grations on 78 sources, with special focus given to covering 
as many of the highest ranked sources as time allowed. Fi- 
nally, in the third campaign, 70 highly ranked sources were ob- 
served with between 60- and 76-minute integrations, with 32 of 
these sources observed multiple times. Observing in this man- 
ner reduced the effects of scintillation and of man-made radio- 
frequency interference (RFI), as most sources were observed 
multiple times with different integration times at different ob- 
serving epochs. 

Data were recorded over 512 filterbank channels of width 
585.3 kHz with a sampling rate of 53 /is[^| Initially all data were 
sampled at 32 bits by the digitisers before being brought down to 
8 bits and written to high-capacity magnetic tape for transporta- 
tion and storage. 

To meet the processing demands imposed by the 5.5 TB of 
data created in the survey, a 22-node 168-core computing cluster 
situated in the MPIfR was used for all data analysis. The PRESTO 
software package (Ransom 2001 1 was used for data processing. 

In the first stage of the processing pipeline, data underwent 
RFI treatment in which a time- and frequency-dependent mask 
was created to be applied at a later stage. As the 1FGL catalogue 
contains no distance information, each beam was de-dispersed 
at 2760 trial dispersion measures in the range 0-997 pc cm -3 
to mitigate against the frequency-dependent delay in the pulsar 
signal due to dispersion by free electrons along the line of sight. 

1 An in-depth review of the observational set-up for pulsar searches 
with the Effelsberg telescope will be published in the initial paper of 
the High Time Resolution Universe North pulsar survey (Barr et al., in 
prep.). 



Our choice of such a fine sampling in dispersion space allows 
for the retention of the data's maximum possible time resolution 
at all dispersion measures. The effect of this was a heightened 
sensitivity to millisecond and potential sub-millisecond pulsars. 

All de-dispersed time series were fast-Fourier-transformed 
(FFT) and the resulting power spectra were de-reddened and 
known RFI frequencies were removed. To reconstruct power dis- 
tributed through harmonics in the Fourier domain, the process 
of incoherent harmonic summing was used. Here the original 
spectra are summed with versions of themselves that have been 
stretched by a factor of two such that all second-order harmonics 
are added to their corresponding fundamental. This process was 
repeated four times such that all power distributed in even har- 
monics up to the 16 th harmonic may be incoherently added to 
the fundamental (see e.g. |Lorimer & K ramer 2005 ). The spec- 
tra from each stage of the summing process were searched for 
accelerated and non-accelerated signals. 

At this stage the number and size of the FFTs required to 
achieve sensitivity to fast binaries becomes computationally too 
expensive for the longest pointings. To deal with this, the pro- 
cessing of the data was split into two stages. Initially all data 
were analysed at full length with a moderate acceleration search 
in the Fourier domain. This analysis is very sensitive to isolated 
and mildly accelerated pulsars in the data. The second process- 
ing stage involves splitting the data into 10-minute blocks and 
re-analysing with a much more intensive acceleration search. Al- 
though this stage uses shorter integrations, it is more sensitive to 
highly accelerated binary systems in the data. The details of the 
acceleration search can be found in Section 4. 

Upon completion of the processing, all candidate signals 
underwent a sifting routine that removes any signals which are 
likely to be RFI. Finally for the top 50 candidates a set of di- 
agnostic plots were created for visual inspection (e.g. Eatough 
et al. 20 10). In cases where there were more than 50 candidates 
with greater than 8-cr significance, the pointing was considered 
to be contaminated by RFI and was flagged for re-observation. 



4 SENSITIVITY 

To estimate the sensitivity of our survey setup we used the ra- 
diometer equation (see e.g. |Lorimer & Kramer 2005 ), 

g . _ p S/Nmi n T ByB ( -Pcyclc \ 5 ,j» 
GyJn p t ohs 8f V 1 ~ Pyclc / 

where the constant factor j3 denotes signal degradation due to 
digitisation, which for 8-bit digitisation is ~ 1%, giving /3 = 
1.01 {Kouwenhoven & Voute 20011. T sys is the system tempera- 
ture of the receiver. From flux density calibration measurements 
we found T sys = 25 K. G is the antenna gain (1.5 KJy -1 at 1.36 
GHz), P C ycio is the pulse duty cycle, t b s is the length of the 
observation, 5f is the effective bandwidth of the receiver (240 
MHz) and n p is the number of polarisations summed, which 
for this survey is always 2. The factor S/N m i„ is the minimum 
signal-to-noise ratio with which we can make a detection. Based 
on false alarm statistics, we chose S/N m i n = 8. 

Figure [T] shows the minimum detectable flux density as a 
function of pulsar period for four dispersion measures; 0, 100, 
500 and 1000 pc cm -3 . Assuming a pulsar of typical duty cycle 
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5 SIMULATIONS 
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Figure 1. Theoretical minimum detectable flux densities (S m i n ) vs. spin 
period for three integration regimes. A minimum detectable signal-to- 
noise of 8 was imposed. The break point at period 10 ms occurs due 
to the assumption that the pulse duty-cycle scales as period -0 5 with a 
maximum value of 1 /3 {Kramer et al.|1998} . 



5% we achieve a minimum detectable flux density of 0.02 mjy 
for a 76-minute pointing and 0.06 mJy for a 10-minute pointing. 

In a binary system, due to the Doppler effect, the appar- 
ent spin frequency of the pulsar drifts with time, spreading the 
pulsar's power in the Fourier domain. To reconstruct Fourier 
power smeared across multiple bins, we employed presto's 
accelsearch routine which uses the 'correlation technique' of 
template matching in the Fourier domain as outlined in |Ran-| 
som et al. (20021. The relationship between the largest range 
of Fourier bins over which to search for drifted signals, Z max , 
and the acceleration to which the search is sensitive, do, can be 
described by 



ao = 



<Pc 



t 2 



(2) 



where P is the spin period of the pulsar being searched for, t b s 
is the integration time and c is the speed of light in a vacuum. 

By assuming a 1-ms period pulsar in a binary system we 
may determine the sensitivity of the search at the extreme limit 
of the pulsar population. For the first processing pass, where 
the data were analysed in full length, we used Z max = 50 to 
achieve sensitivity to accelerations of up to |an| ~42ms" 2 and 
[do | ~ 1 m s~ 2 for 10- and 76-minute integrations respectively. 
For the second processing pass, where data were analysed in 10- 
minute blocks, we used Z max = 600 to achieve sensitivity to 
accelerations of up to |an| ~ 500 m s~ 2 . For all systems which 
fall inside the acceleration limits for their integration lengths, the 
minimum detectable flux density may be calculated via Equation 
□ 



As with any survey that contains many pointings, we may as- 
sume that the chance detection probability for pulsars unassoci- 
ated with the LAT sources targeted is non-negligible. In order to 
determine this probability, simulations of the normal pulsar and 
MSP populations we re made based on the model presented in 



Lorimer et al. 



1 2006 1 using the PSRPOP^Jsoftware. As the popu- 



lation distributions for MSPs and normal pulsars differ, separate 
simulations were performed for pulsars with rotational periods 
above and below 40 ms. 

To simulate the normal pulsar distribution, input model pa- 
rameters were chosen as follows: 

• An empirical period distribution taken from the probability 
density function of the known population. 

• A log-normal luminosity distribution, with mean and vari- 
ance in log space of — 1 . 1 and 0.9, respectively (Faucher-Giguere 
|& Kaspi|2006"l >. 

• A Gaussian distribution of spectral indices, with mean of 
— 1.6 and variance of 0.5. 

• An exponential distribution for the height above the Galac- 
tic plane, with a scale height of 0.33 kpc (|Lorimer et al.|2006) . 

• A radial distribution as described in |Lorimer et al.| (|2006 ). 

• The NE2001 Galactic free electron density model (Cordes 
|& Lazio|200"2) l. 

To simulate the MSP population, the same model param- 
eters were used with the Galactic scale height increased to 0.5 
kpc to better match the known MSP distribution ( |Lorimer et al.| 
[2006) . 

The number of pulsars simulated was such that the cumu- 
lative number of pulsars discovered in simulated versions of the 
Parkes Multibeam Pulsar Survey ( [Manchester et al.|2001 1, the 
Swinburne Intermediate Latitude Pulsar Survey {Edwards et ah] 
2001 1 and its extension ( Jacoby et al. 2009), and the Parkes High 
Latitude Survey ( Burgay et al. 2006 ), matched the real discovery 
numbers. 

Each observation of the 289 sources in the survey was then 
compared against the simulated pulsar distribution to determine 
if a chance detection could have been made. The Galactic pulsar 
population was simulated 1000 times and comparisons repeated 
with each simulation. We find a mean detection rate of 0.4 nor- 
mal pulsars and 0.04 MSPs in the 289 pointings. While not rul- 
ing out the possibility of a chance MSP detection, these figures 
suggest that any MSP found in this survey will be associated to 
the LAT source targeted at 96% confidence. 



6 RESULTS 

6.1 PSR J1745+1017 

6.1.1 Radio analysis 

The main result of this work is the discovery of the radio 
pulsar PSR J1745+1017 in LAT source 1FGL J1745.5+1018 
(2FGL J1745. 6+1015). Initially discovered as a mildly acceler- 
ated candidate in a 10-minute data segment, PSR J 1745+ 10 17 



http://psrpop.phys.wvu.edu/index.php 
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has a period of 2.65 ms and a period derivative of 2.72 x 10~ 21 , 
making it the first MSP discovered with the Effelsberg telescope. 

Upon discovery, PSR J1745+1017 was the subject of an 
intensive timing campaign involving the Effelsberg, Lovell and 
Nancay radio telescopes. Pulse times of arrival (TOAs) were 
analysed with the TEMP02 software package ( Hob bs et al.|2006} 
to create a phase-connected timing solution for the first two 
years of radio observations. The pulsar ephemeris and a selec- 
tion of derived properties are displayed in Table ^ 

From period and period derivative measurements we infer 
a characteristic age of 16.9 Gyr. Combined with the rapid spin 
period of PSR J1745+1017, this implies that the pulsar has un- 
dergone an extensive period of accretion-induced spin-up due 
to mass transfer from a companion star. This hypothesis is sup- 
ported by the timing solution, which shows PSR J 1745+ 10 17 to 
be in a low-eccentricity binary system with an orbital period of 
17.5 hr and a low-mass companion. Using the definition for the 
binary mass function, 

- 2 (asinz) 3 (m,2sini) 3 



f(m 1 ,m 2 ) 



4tt z 



(3) 



G P 2 (mi + m 2 ) 2 ' 

where mi and m,2 are the masses of the pulsar and compan- 
ion respectively; G is Newton's gravitational constant; asini 
is the projected semi-major axis of the system, where i is the 
inclination angle of the binary (defined such that i = 90° is 
edge on); and Pb is the orbital period, we obtain /(mi , TO2) = 
1.4 x 10 _6 Mq. Assuming an average pulsar mass of 1.35 Mq, 
this gives minimum and median companion masses of 0.014 and 
0.016 Mq respectively. 

This small range of low companion masses strongly sug- 
gests that PSR J1745+1017 is a Black Widow system with a 
white dwarf companion that has been heavily ablated by the 
strong particle wind from the pulsar. Unlike Black Widow sys- 
tems such as PSR B 1957+21 ( |Fruchter et aLp988l > and PSR 
J2051-0827 ( |Stappers et al.|1996> , PSR J1745+1017 exhibits 
neither eclipsing behaviour nor DM variations across orbital 
phase, implying a low inclination angle. If the inclination is low, 
then the companion mass is likely to be larger than the presented 
value. As noted by Lorimer & Kramer ( 2005 1, the probability of 
observing a binary system with an inclination less than io for a 
random distribution of orbital inclinations is p(i) = 1 — cos(io). 
This relation suggests that the inclination angle is greater than 
26° at 90% confidence. We therefore place an upper limit of 
0.032 Mq on the mass of the companion to the same confidence 
level. 

From proper motion measurements of PSR J1754+1017, 
we find a transverse velocity of ~ 48 + 9 km s _1 . The small 
dispersion-measure-inferred distance to this pulsar, d ~ 1.3 kpc, 
suggests that the Shklovskii effect (Shklovskii 1970) contribu- 
tion to the measured period derivative will be non-negligible. 
For a pulsar with transverse velocity Vt, the Shklovskii effect 
acts to increase the intrinsic period derivative of the pulsar by a 
factor Pv t jcd (see Cami lo et al.|1994} , where c is the speed of 
light in a vacuum. We find a contribution to the period derivative 
of 5 ± 2 x 10 -23 or about 18% of the measured value. 

Assuming a moment of inertia of 10 45 g cm 2 and using 
the Shklovskii-corrected period derivative, we derive the spin- 
down luminosity of PSR J 1745+ 1017 to be E = 4ty 2 IP/P 3 = 
4.7 x 10 33 ergs s _1 , a value which makes this pulsar a good 
candidate for pulsed 7-ray emission (Abdo et al.|2010a| >. 



Table 1. PSR J1745+1017 ephemeris created from TOAs taken with the 
Nancay, Effelsberg and Lovell telescopes over the course of 22 months. 
Numbers in parentheses represent twice the formal l-cr uncertainties in 
the trailing digit as determined by TEMP02. The dispersion-measure- 
derived distance was estimated using the NE2001 Galactic electron den- 
sity model (Cordes & Lazio 2002J. Due to intrinsic uncertainties in this 
model, the estimation is likely to have an uncertainty of ~ 20%. The 
mass function calculation assumes an average mass of 1.35 Mq for the 
pulsar. The characteristic age, spin-down luminosity and surface mag- 
netic field strengths were calculated using the Shklovskii-corrected pe- 
riod derivative. The position, frequency and DM are all measured with 
respect to the given reference epoch. These parameters were determined 
with TEMP02, which uses the International Celestial Reference System 
and Barycentric Coordinate Time. Refer to Hobbs et al. ( 2006 1 for infor- 
mation on modifying this timing model for observing systems that use 
TEMPO format parameters. 

PSR J1745+1017 ephemeris 



Fitted timing parameters 



- 1 ) 



Right Ascension (RA. J2000) (hh:mm:ss) 
Declination (Decl. J2000) ( ° :':") 
Proper motion: 

in RA. (fi a cos(Decl)) (mas yr~ 
in Decl. (fig) (mas yr -1 ) 
Period (s) 

Period derivative (x 10 -21 ) 
Dispersion measure (pc cm -3 ) 
Orbital period (days) 
Projected semi-major axis (lt-s) 
Epoch of ascending node (MJD) 
re (= e cosw) (x 10~ 5 ) 
T] (= esinw) (xl0~ 5 ) 

Fixed parameters 



17:45:33.8371(7) 
+10:17:52.523(2) 

6(1) 
-5(1) 

0.00265212967108(3) 
2.73(1) 
23.970(2) 
0.730241444(1) 
0.088172(1) 
55209.968794(2) 
0(2) 
0(2) 



Reference epoch (MJD) 
Clock correction procedure 
Time system 

Solar system ephemeris model 
Binary model 

Derived parameters 



55400 
TT(TAI) 
TCB 
DE414 



ELL1 (Lange et al. 



2001 



Frequency (Hz) 

Frequency derivative (Hz s _1 X 10 -16 ) 

Orbital eccentricity (x 10~ 5 ) 

Epoch of periastron passage (MJD) 

Galactic longitude (J2000) (deg) 

Galactic latitude (J2000) (deg) 

Mass function (Mq ) 

Minimum companion mass (Mq ) 

Median companion mass (Mq ) 

Dispersion measure-derived distance (kpc) 

Shklovskii-corrected 



period derivative (x 10 



Characteristic age (Gyr) 
Spin-down luminosity (x 10 33 erg s — x ) 
Surface magnetic field (x 10 7 G) 
rms residual (ps) 

Further parameters 



377.05547013813(5) 
-3.88(2) 
0(2) 
55210.3(4) 
34.8693081(3) 
19.2536887(5) 
1.38(l)xl0- 6 
0.0137 
0.0158 
1.3(2) 

2.22(5) 
18.9 
4.7 
7.7 
5.05 



Median flux density at 1.36 GHz (mJy) 
Maximum flux density at 1.36 GHz (mJy) 
Span of timing data (MJD) 
Number of TOAs 



55225 



0.3 
4.4 
■ 56026 
156 
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6.1.2 Gamma-ray analysis 

In order to characterize the 7-ray emission of PSR J 1745+1017, 
we selected Fermi LAT data recorded between 2008 August 4 
and 2012 March 7, with reconstructed energies larger than 0.1 
GeV, directions within a circular region of interest (ROI) of 15° 
radius around the pulsar's position, and zenith angles smaller 
than 100°. We further restricted the dataset to "Source" class 
events of the P7_V6 instrument response functions, and rejected 
times when the rocking angle of the LAT exceeded 52° or when 
the Earth's limb infringed upon the ROI. The selected 7-ray 
events were finally phase-folded using the ephemeris given in 
Table [T] and the Fermi plug-in distributed with TEMP02 ( |Ray| 
|etal.|2011| >. 

Initial pulsation searches using standard data selection cuts 
yielded marginal detections only. For instance, selecting photons 
found within 1 of the pulsar and with energies larger than 0. 1 
GeV, we found an //-test parameter ( de Jager & Buschi ng|2010[ > 
of 15.6, which translates to a significance of ~ 3.1a. Nonethe- 
less, pulsation searches can be made more sensitive by weight- 
ing the photons by the probability that they originate from the 
pulsar. These probabilities can be computed through a spectral 
analysis of the pulsar and the neighbouring sources l |Kerr|2011| 
IGuillemot et al.|2012| >. 

The 7-ray spectrum of PSR J1745+1017 was measured by 
fitting sources in the ROI using a binned likelihood method, 
with the pyLikelihood module included in the Fermi Science 
Tool0 The source model used for the analysis included the 
spectral parameters of the 78 sources of the Fermi LAT Sec- 
ond Source Catalogue (2FGL Nolan et al. 20121 found within 
20° of the pulsar. The spectrum of PSR J1745+1017 was rep- 
resented by an exponentially-cut-off power-law of the form 
dN/dE = N (E/l GeV)" r exp(-_B/£ c ), where N is a 
normalization factor, F is the photon index, and E c is the cut- 
off energy. The extragalactic diffuse emission and the residual 
instrument background were modelled using the iso_p7v6source 
template, and the Galactic diffuse emission was modelled us- 
ing the gal^2yearp7v6jv0 map cube. The spectral parameters 
of the sources within 5° of the pulsar as well as the normal- 
izations of the diffuse models were re-fit, while the parameters 
of other sources were fixed at the 2FGL catalogue values. We 
found no evidence of significant emission from the pulsar in the 
phase range [0.6; 1]. In order to increase the signal-to-noise ra- 
tio of the pulsar we thus restricted the dataset to photons with 
reconstructed pulse phases between and 0.6 (the 7-ray light 
curve of PSR J1745+1017, presented below, is indeed compati- 
ble with showing emission only in this phase range). The best- 
fitting spectral parameters of PSR J1745+1017 are listed in Table 
[2] along with the integrated photon and energy fluxes above 0. 1 
GeV derived from these results. Systematic uncertainties were 
calculated by running the same analysis as described above, but 
using bracketing IRFs for which the effective area has been per- 
turbed by ±10% at 0.1 GeV, ±5% near 0.5 GeV, and ±10% at 
10 GeV, with linear interpolations in log space between. 

Using the best-fitting spectral model for the ROI obtained 
from the analysis described above, we could calculate probabili- 
ties that the photons in the ROI originate from PSR J1745+1017. 

3 http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.html 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
Pulse Phase 



Figure 2. Multi-wavelength light curves of PSR J1745+1017. The bot- 
tom panel shows the 1.4 GHz radio profile recorded with the Effels- 
berg Radio Telescope. The upper panels show probability-weighted 7- 
ray light curves in different energy ranges, with 30 bins per rotation. 
Horizontal dashed lines show the estimated 7-ray background levels. 
Two rotations are shown for clarity. 



Selecting events found within 5° and with calculated probabili- 
ties larger than 0.01, we obtained a weighted //-test parameter 
(see Kerr 2011) of 58.0, corresponding to a pulsation signifi- 
cance of 6.5a. Figure[2]shows probability-weighted light curves 
for PSR J1745+1017 in different energy bands. The background 
levels shown in Figure [2] were calculated by summing the prob- 
abilities that selected 7-ray events do not originate from the pul- 
sar, as described in Guillemot et al. (2012). Statistical error bars 
were obtained by calculating y/2~2i w h where Wi represents the 
photon probability and i runs over photons falling in the same 
phase bin ( [Pletsch et al.|2 012 1. As can be seen from Figure|2] the 
7-ray profile of PSR J1745+1017 shows evidence for two dis- 
tinct peaks, at phases ~ 0.14 and ~ 0.39. Fits of the two 7-ray 
peaks with Lorentzian pulse shapes above 0.3 GeV yielded the 
peak positions <I>; and the full widths at half maxima FWHJVL 
listed in Table[2] We also attempted to fit the light curve with one 
asymmetric Lorentzian pulse shape, and found that the model 
with two peaks is slightly preferred, at the ~ lcr level. The radio- 
to-7-ray lag 5 = $1 — $ r (where <J> r = is the phase of the 
maximum of the radio profile shown in Figure [2} and the 7-ray 
peaks separation A = $2 — $1 are found to be S = 0.14 ±0.04 
and A = 0.26±0.06 (see Table[2]l. The uncertainty on the radio- 
to-7-ray lag due to the error on the measurement of the DM 
parameter is estimated to be A(5) = A(DM) / (K f 2 ) where 
K = 2.410 x 10~ 4 MHz -2 cm" 3 pc s -1 is the dispersion 
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Parameter 



Value 



First peak position, $1 0.14 ±0.04 

First peak full width at half-maximum, FWHMi 0.20 ±0.16 

Second peak position, $2 0.39 ± 0.03 

Second peak full width at half-maximum, FWHM 2 .... 0. 10 ± 0. 10 

Radio-to-7-ray lag, 8 0.14 ±0.04 

7-ray peak separation, A 0.26 ± 0.06 

Photon index, T 1.6 : 0.2 '^j 

Cutoff energy, E c (GeV) 3.2 ± 1.2 1 

Photon flux (> 0.1 GeV), F100 (10~ 8 cm" 2 s" 1 ) .... 1.1 ± 0.3to's 

Energy flux (> 0.1 GeV), G100 (10~ 12 erg cm" 2 s" 1 ) 9.3 ± 1.2±g;| 

Luminosity, L 7 / f n (10 33 erg s" 1 ) 1.8 ± 0.6±o.e 

Efficiency, 77 / f n 0.3 ± 0.1±%{ 



Table 2. Measured 7-ray light curve and spectral parameters for PSR 
11745+1017. First quoted uncertainties are statistical, and the second are 
systematic. Details on the determination of the systematic uncertainties 
are given in Section [6. 1 ,2| 



constant. We find A(<5) ~ 10~ 3 x P, which is very small com- 
pared to the statistical error bar. Such values of 8 and A are rel- 
atively common amongst other known 7-ray pulsars (see Figure 
4 of JAbdo et al.|2010b) , and match the predictions of theoretical 
models that place the high-energy emission from pulsars at high 
altitudes in the magnetosphere ( Romani & Yadigaroglu ) 1995| >. 

The ephemeris used for phase-folding the 7-ray data con- 
sidered in this analysis is based on radio timing taken after 2010 
January 30 (MJD 55226). In an attempt to determine whether the 
ephemeris describes the rotational behaviour of the pulsar over 
the entire Fermi LAT dataset accurately, we analysed the evo- 
lution of the weighted //-test parameter as a function of time. 
The results of this analysis are shown in Figure [3] First of all, 
because the //-test depends linearly on the number of photons 
for a given pulsed signal fraction, the linear increase of the li- 
test parameter as a function of the dataset length provides further 
evidence that the pulsed 7-ray signal from PSR J1745+1017 is 
real. Moreover, the increase of the //-test parameter when going 
forward or backward in time is monotonic outside the formal 
ephemeris validity interval, which indicates that the ephemeris 
given in Table [T] provides a good description of the pulsar's ro- 
tational behaviour across the entire LAT dataset. 

From the energy flux G100 measured from the spectral 
analysis, we calculated the 7-ray luminosity above 0. 1 GeV us- 
ing L-y = 4irfnGiood 2 , where fn is a geometrical correction 
factor depending on the beaming angle of the pulsar and the 
viewing geometry (Wa tters et al.|200"9) . Using the distance de- 
rived from the NE2001 model of d = 1.3±0.2 kpc, we get the 
7-ray luminosity L 7 //n ~ 1.8 x 10 33 erg s _1 and the 7-ray ef- 
ficiency, r\l fn — L-y/E ~ 0.3. This value for the efficiency lies 
well within the distribution of 7-ray efficiencies seen for other 
LAT-detected MSPs ( |Abdo et al.|2010bl >. 

From the measurement of the 7-ray energy flux and the 
proper motion of the pulsar, we can put an upper limit on the 
distance derived from the inequality L 7 < E (see Equation 2 of 
Guillemot etal.|2012|. Assuming / = 10 45 g cm 2 and fn = 1 (a 




54600 54800 55000 55200 55400 55600 55800 56000 

Time (MJD) 



Figure 3. Evolution of the weighted //-test Test Statistic (TS) as a func- 
tion of time. The blue, dashed line shows the weighted //-test parameter 
with increasing time, using data taken from 2008 August 8. The red, dot- 
ted line shows the weighted //-test parameter when going backwards in 
time, using data taken before 2012 March 7. The vertical, dashed line 
indicates the formal start of the ephemeris validity interval. 



6.2 Radio pulsar non-detections 

During the course of the observations, other observatories dis- 
covered six radio-emitting pulsars associated with the sources 
observed in this survey. These pulsars were not detected in this 
survey due to a combination of reasons, such as limited flux den- 
sity, local RFI conditions and source position redefinitions. Fur- 
ther information on all sources observed in this work may be 
found in the corresponding on-line material, with the 15 highest 
ranked sources presented in Table |Al| 



6.2.1 PSRJ2030+3641 



typical value for other 7-ray MSPs, see e.g. jVenter et al.|2 009 1, 
we get d 1.9 kpc. 



PSR J2030+3641 < |Camilo et ^[2012) , initially undetected in 
our data processing, exhibits a rotational frequency which is 
similar to the 10 th sub-harmonic of the 50-Hz local mains fre- 
quency. Unfortunately signals which are similar in frequency to 
harmonics of the mains are often lost in noise or excised from 
the data in the process of Fourier RFI excision, wherein known 
RFI signals are removed from the power spectrum of the obser- 
vation prior to pulse detection. It should be noted that this pulsar 
was discovered using the Robert C. Byrd Green Bank Telescope 
in the US, where the mains frequency is 60 Hz. 

Reprocessing of these data with the published ephemeris 
yielded a detection of the pulsar at signal-to-noise 6.4, a value 
which lies below the detection threshold of this survey. Using 
Equation[T] we found a corresponding estimated flux density of 
0.08 mjy. This value is half that of the published flux density 
at 1.5 GHz, a discrepancy which may be attributed to the ex- 
tra noise induced by the oscillations of the mains electricity at 
the folding period. It should be noted, that the derived position 
of J2030+3641 is coincident with the position observed in this 
survey. 
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PSR 


Pointing 


•Srriin 


Reference 




offset (°) 


(mJy) 




J0734-1559 


0.05 


0.1 


Saz Parkinson 1 201 1 




J1803-2I49 


0.05 


0.1 


|Pletsch efaqpHI 




J2030+4415 


0.02 


0.08 


|pietsch efaT|pjn 




J2 139+47 16 


0.04 


0.09 


|pietsch eTaT|pjn 





Table 3. Limiting flux densities (.S m j n ) for LAT detected pulsars coin- 
cident with sources observed in this work. All pulsars are assumed to 
have a radio pulse with a 10% duty-cycle and scattering effects are ig- 
nored. It should be noted that the upper radio flux limits presented here 
are higher than those published in Pletsch et al. ( 20 1 2 1, due to the use of 
an increased detection threshold of signal-to-noise 8. 



6.2.2 PSRs J1646-2142, J 1816+4510 and J 1858-2218 



PSRs J1646-2142, J1858— 2218 ( |Ray et al.| [20T2> and 
Jl 8 16+45 10 (Kaplan et aLl[2012) were all discovered at LAT 
positions that were significantly different from the position de- 
termined in the 1FGL catalogue. Because of this, the discovery 
position of these pulsars is a half-beamwidth or more away from 
the position observed in this work, and so no detection is ex- 
pected. 

6.2.3 PSRs J0307+7443 and J1828+0625 



PSRs J0307+7443 and J1828+0625 ( |Ray et al.|2012) were both 
discovered at positions coincident with sources observed in this 
survey. Both pulsars were detected at low radio frequency and 
appear to be weak radio emitters. Assuming similar profile char- 
acteristics at both 1.36 GHz and at the discovery frequency, we 
found upper limits of 0. 1 and 0.2 mJy on the radio flux densities 
at 1.36 GHz for PSRs J0307+7443 and J1828+0625, respec- 
tively. 

6.3 Gamma-ray pulsar non-detections 

Through blind searches of Fermi LAT photons, a further four 
pulsars associated with the sources observed have since been 
discovered (Pletsch et al.|2012|rSaz Parkinson|2011) . No radio 
detection has so far been made for these pulsars, therefore in Ta- 
ble [3] we present our upper limits on their radio flux densities at 
1.36 GHz. 

A simple model for the telescope beam was used to adjust 
limiting flux densities for position offsets between the pointing 
and true position of the source. Limiting flux densities were di- 
vided by a factor q, given by q = e - ' 6 /^ / 15 , where is the 
pointing offset and (f> is the beam half-width at half maximum, 
to return the true limiting flux density at 1.36 GHz (S m i n ). 



7 DISCUSSION 

Considering the success that other surveys of unassociated Fermi 
LAT sources have h ad (Keith et al.|2011| |Cognard et al.|2011| 



Rans om et al.|20lT] |Camilo et al.|2012*| T we must consider the 
reasons behind our relatively small yield of discoveries. 

Sources observed in this work were selected based on infor- 
mation provided in the 1FGL catalogue. However, in the recently 



published 2FGL catalogue, many of these sources are no longer 
listed. A detailed description of the possible causes for this can 
be found in Section 4.2 of |Nolan et"al] ( |2012| >. Using the 2FGL 
catalogue to examine the 289 sources selected for observation, 
we find: 

• 106 sources that are no longer listed. 

• 46 sources that have moved by more than half the FWHM 
of the telescope beam. 

• 22 sources marked as potentially spurious. 

Selecting against these sources and choosing only sources 
for which the entire updated 95% confidence region was con- 
tained within the FWHM of the telescope beam, we are left with 
72 sources for which observations adequately cover the LAT 
source targeted. 

This would suggest that the detection rate for the survey is 
~1.4%. However, as noted by Camilo et al. (2012), a large num- 
ber of pulsar surveys at radio wavelengths have focused on ob- 
serving the Galactic plane at low latitudes, introducing a strong 
bias against the discovery of new radio-emitting pulsars in this 
volume. Selecting only those 58 sources that lie in the region 
\b\ > 3.5°, we find the detection rate to be 1.7%. 

We note that the strategy employed by other surveys of LAT 
unassociated sources has generally been to select a small num- 
ber of sources and observe with long integration times, often 
at lower frequencies (e.g. Ransom et al.|20lT) . By contrast, we 
have selected a large number of sources and observed at various 
integration lengths at higher frequency. Although this strategy 
has lacked the discoveries that other searches have made, it has 
provided the most complete picture yet of the LAT sky at high 
radio frequencies. 

The results of our survey suggest that future surveys of 
unassociated LAT sources would benefit from lower observing 
frequencies, as wider beam-widths will increase the LAT er- 
ror ellipse coverage, which, when combined with the decrease 
in LAT error ellipse size due to increasing observation depth, 
will significantly reduce the problems associated with source re- 
definitions. Furthermore, we note that the majority of pulsars 
discovered through radio searches of the LAT data have been 
found at low frequencies. This suggests that this population of 
pulsars may have particularly steep spectral indices. Contin- 
ued multi-frequency study of pulsars discovered in unassociated 
LAT sources will deny or confirm this hypothesis. 



8 CONCLUSION 

We have performed repeated observations of 289 unassociated 
sources from the First Fermi LAT catalogue of 7-ray sources, 
resulting in the discovery of PSR J1745+1017. 

PSR J1745+1017 has a 2.65-ms spin period and orbits a 
very low-mass companion with an orbital period of 17.5 hours. 
The low mass of the companion suggests that PSR J 1745+ 10 17 
is a Black Widow system, however a two-year timing campaign 
has thus far found no evidence of eclipses, dispersion measure 
variations across orbital phase or Shapiro delay. By performing 
weighted-probability analysis of LAT photons in the region of 
interest around PSR J1745+1017, we also detected 7-ray pul- 
sations from the source. This gives a clear indication that PSR 
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J1745+1017 is responsible for the 7-ray emission seen from 
1FGL J1745.5+1018 (2FGL J1745.6+1015). 

The observed sample contains three radio pulsars, newly 
discovered through similar searches performed by other ob- 
servatories, and four 7-ray selected, radio-quiet pulsars, found 
through blind searches of Fermi LAT data. 

No radio detection was made for the four 7-ray selected 
pulsars contained in our sample. The upper limits on their radio 
flux densities at 1.36 GHz were calculated. These limits can be 
found in Table[3] 

Of the three newly discovered radio pulsars coincident with 
LAT sources searched in this survey, only PSR J2030+3641 was 
detected using the known timing solution. PSRs J0307+7443 
and J 1828+0625 were suspected to have radio flux densities at 
1.36 GHz lower than the limiting flux density of our observa- 
tions. 
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